KANAZAWA 99-26 



An Almost Perfect Quantum Lattice Action for Low-energy 

SU(2) Gluodynamics 



Maxim N. Chernodub a , Shouji Fujimoto 6 , Seikou Kato fe , Michika Murata 6 , 
Mikhail I. Polikarpov a and Tsuneo Suzuki 6 
a ITEP, B.Cheremushkinskaya 25, Moscow 117259, Russia 
b Institute for Theoretical Physics, Kanazawa University, Kanazawa 920-1192, Japan 

(February 1, 2008) 



Abstract 

We study various representations of infrared effective theory of SU (2) Glu- 
odynamics as a (quantum) perfect lattice action. In particular we derive a 
monopole action and a string model of hadrons from SU (2) Gluodynamics. 
These are lattice actions which give almost cut-off independent physical quan- 
tities even on coarse lattices. The monopole action is determined by numerical 
simulations in the infrared region of SU (2) Gluodynamics. The string model 
of hadrons is derived from the monopole action by using BKT transforma- 
tion. We illustrate the method and evaluate physical quantities such as the 
string tension and the mass of the lowest state of the glueball analytically 
using the string model of hadrons. It turns out that the classical results in 
the string model is near to the one in quantum SU(2) Gluodynamics. 

PACS: 12.38.Gc, 11.15.Ha 
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I. INTRODUCTION 



Low-energy effective theory of QCD is important for analytical understanding of hadron 
physics. Before derivation of such an effective theory we have to explain the most important 
non-perturbative phenomenon, quark confinement. Wilson's lattice formulation [|TJ shows 
that the confinement is a property of a non-abelian gauge theory of strong interaction. 
At strong coupling the confinement is proved analytically. At weak coupling (near to the 
continuum limit) there are a lot of numerical calculations showing the confinement of color. 
The mechanism of confinement is, however, still not well understood. One of approaches to 
the confinement problem is to search for relevant dynamical variables and to construct an 
effective theory in terms of these variables. 

From this point of view the idea proposed by 't Hooft || is very promising. It is based 
on the fact that after a partial gauge fixing (abelian projection) SU(N) gauge theory is 
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reduced to an abelian [/(l)^ 1 theory with N — 1 different types of abelian monopoles. 
Then the confinement of quarks can be explained as the dual Meissner effect which is due 
to condensation of these monopoles. The QCD vacuum is dual to the ordinary supercon- 
ductor: the monopoles playing the role of the Cooper pairs. The confinement occurs due 
to the formation of a string with electric flux between quark and anti-quark. It is a dual 
analogue of the Abrikosov string ||. The mechanism of confinement is usually called the 
dual superconductor mechanism. 

There are many ways to perform the abelian projection, but in the Maximal abelian (MA) 
gauge Q many numerical results support the dual superconductor picture of confinement 
in the framework of lattice Gluodynamics (see, for example, reviews |||7J]). These results 
suggest that the abelian monopoles which appear after the abelian projection of QCD, are 
relevant dynamical degrees of freedom in the infrared (IR) region. We expect hence, after 
integrating out all degrees of freedom other than the monopoles, an effective theory described 
by the monopoles works well in the IR region of Gluodynamics. 

The effective monopole action on the MA projection of SU(2) lattice Gluodynamics was 
obtained by Shiba and Suzuki || using an inverse Monte-Carlo method [[J. Assuming 
that the lattice action contains only quadratic terms of monopole currents, they found that 
the action has a form theoretically predicted by Smit and van der Sijs (TIJ. This was the 
first derivation of an effective theory of lattice Gluodynamics in terms of the monopole 
currents. However the steps of block-spin transformation performed in || were rather few 
to see the continuum limit. In |TT| they considered also four- and six-point interactions 
assuming a direction symmetric action on the large (48 4 ) lattice. More steps of the block-spin 
transformations were carried out also. It is stressed that the action seems to satisfy a scaling 
behavior, that is, it depends on the physical length b = na{f3) alone, where n is the number 
of the blocking transformations and a((3) is the lattice spacing. This remarkable scaling is 
consistent with the behavior of the perfect action on the renormalized trajectory (RT) which 
is an effective theory in the continuum limit formulated on the lattice with the lattice distance 
b. Here b plays a role of the physical scale at which the effective theory is considered. On 
RT, although we can predict physical quantities only on the b lattice sites, they are the same 
as evaluated from the continuum theory. For example, the continuum rotational invariance 
should be satisfied. The restoration of the continuum rotational invariance for the quark- 
antiquark static potential was studied using a naive Wilson loop operator. However, the 
continuum rotational invariance was not confirmed in the IR region of SU(2) Gluodynamics 
12] . This is because the cut-off effect of such an operator is of order of the lattice spacing 



of the coarse lattice. To check restoration of the continuum rotational invariance, we should 
determine the correct form of physical operators (the perfect operator) as well as the perfect 
action on the blocked lattice. 

The main task of this publication is to derive the perfect monopole and the string action 
as an low-energy effective theory of SU(2) Gluodynamics and evaluate physical quantities 
analytically using a renormalized operator. 

In Section [TI] we discuss how to derive the renormalized monopole and the string action 
from SU(2) Gluodynamics. We show new results of the analysis of the monopole action 
which is obtained by using inverse Monte-Carlo method. In Section [TT1] we discuss how to 
construct the perfect operator for the static potential. In Section [TV] we calculate the string 
tension and the glueball mass for the SU(2) Gluodynamics in terms of the strong coupling 
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expansion of the string model analytically. It turns out that the classical results in the 
string model is near to the one in quantum SU(2) Gluodynamics. The continuum rotational 
invariance of the static potential is shown also analytically. In Section [V] we analyse the 



numerical results in details. Section [VI] is devoted to concluding remarks. 



II. ALMOST PERFECT MONOPOLE ACTION FROM SU(2) GLUODYNAMICS 

A. Our method 

The method to derive the monopole action is the following: 

1 We generate SU{2) link fields {U(s,/i)} using the simple Wilson action for SU(2) 
Gluodynamics. We consider 24 4 and 48 4 hyper-cubic lattice for (3 = 2.0 ~ 2.8. 

2 Next we perform an abelian projection in the Maximal abelian gauge to separate 
abelian link variables — 7T < 0fj,(s) < tt) from gauge fixed SU(2) 
link fields. 

3 Monopole currents can be defined from abelian plaquette variables 9^ u (s) following 
DeGrand and Toussaint [TBI. The abelian plaquette variables are written by 



^(s) = ^(s) + ^(s + //)-^( S + i>)-^(s), (-4tt < 9^{s) < 4tt). (1) 
It is decomposed into two terms: 

6^(s) = V(s) + 27T7V(S), (-7T < 9^(s) < Tf). (2) 

Here, 0^ u (s) is interpreted as the electro-magnetic flux through the plaquette and the 
integer n^ u (s) corresponds to the number of Dirac string penetrating the plaquette. 
One can define quantized conserved monopole currents 

k^s) = -e^padvUpais + ft), (3) 

where d denotes the forward difference on the lattice. The monopole currents satisfy a 
conservation law d'^k^(s) = by definition, where d' denotes the backward difference 
on the lattice. 

4 We consider a set of independent and local monopole interactions which are summed 
up over the whole lattice. We denote each operator as Si [k] . Then the monopole action 
can be written as a linear combination of these operators: 

S[k] = J2G l S l [k}, (4) 

i 

where Gi are coupling constants. 
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We determine the set of couplings Gi from the monopole current ensemble {k^s)} with 
the aid of an inverse Monte-Carlo method first developed by Swendsen and extended 
to closed monopole currents by Shiba and Suzuki [|8]|J. 

Practically, we have to restrict the number of interaction terms. It is natural to assume 
that monopoles which are far apart do not interact strongly and to consider only short- 
ranged interactions of monopoles. The form of actions adopted here is 27 quadratic 
interactions and 4-point and 6-point interactions. We have not assumed a direction 



symmetric form of the action as done in The detailed form of interactions are 



shown in Appendix A. Note that all possible types of interactions are not independent 
due to the conservation law of the monopole current. So we get rid of almost all the 
perpendicular interactions by the use of the conservation rule. The validity of the 
truncation has been studied and supported in the earlier works. For details, see []8,iT 



5 We perform a block-spin transformation in terms of the monopole currents on the 
dual lattice to investigate the renormalization flow in the IR region. We adopt n = 
1, 2, 3, 4, 6, 8 extended conserved monopole currents as an n blocked operator ||14|| : 



n-l 



KM n) )= E k^ns^ + (n-l)fi + tO + jp + la) (5) 

i,j,l=0 

= %(* (n) ), (6) 

The renormalized lattice spacing is b = na(f3) and the continuum limit is taken as the 
limit n — > oo for a fixed physical length b. 

We determine the effective monopole action from the blocked monopole current ensem- 
ble Then one can obtain the renormalization flow in the coupling constant 
space. 

The physical length b = na(/3) is taken in unit of the physical string tension y/<7 phys . 
We evaluate the string tension <jL at from the monopole part of the abelian Wilson 
loops for each f3 since the error bars are small in this case. The lattice spacing a{(3) is 
given by the relation a(/3) = J 1 aiatl 'vphys [0]. Note that b = I.Oct"^ 2 corresponds to 
0.45/m, when we assume a phys = (440MeV r ) 2 . 



B. Numerical results 



We list new results below in comparison with earlier numerical analysis of the monopole 
action. 

1 The inverse Monte-Carlo method works well and the coupling constants of the action 
are fixed beautifully. The quadratic coupling constants and 4-point coupling con- 
stant are plotted versus the physical length b = na(f3) for each n extended monopole 
in Figure [3]. The first three figures show quadratic self coupling Gi(b), quadratic 
nearest-neighbor couplings (G 2 (b) (black symbol), G 3 (b) (open symbol)) and G w (b), 
respectively. The self-coupling term is dominant and the coupling constants decrease 
rapidly as the distance between the two monopole currents increases. 
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Gi(6) » G 2 (b) ~ G 3 (6) > • • • > Gio(6) > • • • 



The 4-point coupling constant becomes negligibly small in comparison with the 

1 Icy 

quadratic couplings for large b region (6 > l.ba ph y S ). The 6-point coupling constant 
behaves similarly as the 4-point coupling does and becomes much smaller for large b 
region. 

quadratic couplings >> 4-point coupling >> 6-point coupling 



^From these figures we see a scaling of the action iS[/c M , n, a(/3)] — > S[K^, b = na(/3)} 
for fixed physical length b = na{(3) looks almost good for n > 4. The obtained action 
appears to be a good approximation of the action on the RT. 

In Figure we plot the projected lines (Gi(b) — G 2 (b), G 2 (b) — G 3 (b) and G\(b)-A- 
point, respectively) of the renormalization flow. Each flow line for smaller (3 ( which 
corresponds to larger b ) is beautifully straight with very small errors. The quadratic 
interactions for monopoles are dominant for larger b , that is, only the quadratic 
interaction subspace seems sufficient in the coupling space for low-energy SU(2) Gluo- 
dynamics. We also see the effective monopole action tends to go to the weak coupling 
region when we go to the infrared region of SU(2) Gluodynamics. 

The quadratic coupling constants at b = 2.14 are plotted versus the squared distance 
R 2 in unit of squared physical length b 2 in Figure |3|. We see the direction asymmetry 
of the current action. (For example, G 2 ^ G3.) This behavior of the action does not 
occur in the case of compact QED, because the monopole action can be obtained from 
the Villain form of compact QED exactly in an analytical way and it does not depend 
on the direction between two monopole currents. In jlT) they have neglected this effect 



and have considered a direction symmetric form of the monopole action but as we will 
see later that this direction asymmetry of the current action is natural and important 
features of the perfect lattice action. 



III. A PERFECT OPERATOR FOR PHYSICAL QUANTITIES 



In previous section we have studied the renormalized monopole action S[k] performing 
block spin transformation up to n = 8 numerically, and have found the scaling for fixed 
physical length b looks almost good. If the continuum rotational invariance of physical 
observables is satisfied in addition in the framework of S[k], we can regard S[k] as a good 
approximation of RT. 



A. Improved and perfect operator 



In Gluodynamics, the string tension from the static potential is one of important physical 
quantities. However, it is a problem how to evaluate the static potential between electrically 
charged particles after abelian projection. In the earlier work we considered a naive 
abelian Wilson loop operator and S[k] on the coarse lattice to evaluate the static potential, 
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but the continuum rotational invariance of the potential could not be well reproduced even 
for the infrared region of SU(2) Gluodynamics. This is because the cut-off effect of such an 
operator is of order of the lattice spacing of the coarse lattice. Only the scaling behavior of 
the action is insufficient. We should also adopt improved physical operators on the coarse 
lattice in order to get the correct values of physical observables. An operator giving a cut-off 
independent value on RT is called perfect operator. 



B. The method 



As will be shown in Subsection III D , when we consider a monopole action composed of 



general quadratic interactions alone, a block spin transformation can be done analytically 



T5| . We find a perfect operator for a static potential starting from an operator in the 
continuum limit. The continuum rotational invariance is shown exactly with the operator. 
This is an example of a perfect operator. 

What happens in low-energy SU(2) Gluodynamics? It is natural that one can not per- 
form a block spin transformation analytically. However, as shown in the previous section, 
the abelian monopole action S[k] which is obtained numerically is well approximated by 
quadratic interactions alone for large b. The monopole action on the renormalized trajec- 
tory (RT) is expected to be near to the quadratic coupling constant plane in the infrared 
region. We can perform the analytic block spin transformation along the flow projected on 
the quadratic coupling constant plane as shown in Figure [|. When we define an operator on 
the fine a lattice, we can find a perfect operator along the projected flow in the a — > limit 
for fixed b. Let us adopt the perfect operator on the projected space as an approximation of 
the correct operator for the action S[k] on the coarse b lattice. It will be shown in the fol- 
lowing Subsection [TirH that the above standpoint may be justified as long as the quadratic 
monopole interactions are dominant. 



C. Various operators for a static potential 

There is another problem what is the correct operator for the abelian static potential 
in abelian projected SU(2) Gluodynamics on the fine a lattice. First let us consider the 
following abelian gauge theory of the generalized Villain form on a fine lattice with a very 
small lattice distance: 

S[e,n] = j~2 E (d [ti e u] (s) + 2Tm ia ,(s)){A L D )(s- s')(d [fl 9 v] (s') + 27rn^(s')), (7) 

where # M (s) is a compact abelian gauge field and the integer- valued tensor n^ v {s) comes from 
the periodicity of the lattice action ([/]). Both of the variables are defined on the original 
lattice. Al(s-s') = —dd 5 S)S > is the lattice Laplacian and we write D = (3Aj} + D for later 
convenience, where D is a general operator. Since we are considering a fine lattice near to 
the continuum limit, we assume the direction symmetry of D . Note that D = 2n 2 j3yA~[} 
corresponds to the ordinary Villain action for compact QED. In this type of model, it is 
natural to use an abelian Wilson loop W{C) = expi X)e(^( s )> ^( s )) f° r particles with 
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fundamental abelian charge, where J M (s) is an abelian integer-charged electric current. The 
expectation value of W(C) is written as 

(W(C)) = /exp(z£ J,(s)e^s)\) = Z[J]/Z[0], (8) 



Z[J]= UdB^s) Yl ex P \-s[e,n]+iJ2J,(s)W\ 

n n uv (s)=-oo I J 



(9) 



*pi/(f)= 

Next it is known that the theory with the above action (]?]) is equivalent to the lattice 



form of the modified London limit of the dual abelian Higgs model |T6| as shown in Appendix 
B 

«S[C,«M = i- £ {d { ,C v] {s)f 

^P s;pL>v 

+ \ E (Wis) - C,(s) + 2irl,(s))D' - 1 (s - s')(<W) - C,(s') + 2^{s')). (10) 

s,s';fi 

The static potential for electrically charged particles is evaluated by a dual 't Hooft operator 
H(C) = exp j-i- £ {d^C v] {8) - 2^S-l u {s)f + ^ E (^](s)) 2 } , (H) 

I ^P s;fi>u ^P s;/i>v 

where *S^(s) is dual to the surface which is spanned inside the contour J^s). 



Thirdly, when use is made of the BKT transformation [|i7| -|l"9"lj , the action (0) is equivalent 
to the following monopole action 

S[K{s)\ = E Us)D (s - s')k^s'). (12) 

s,s' ,fl 

We see that the area law term is given correctly also by the following operator in the 
monopole representation as shown in Appendix B: 

W m (C) = exp ^TrzX^OOA^s)), (13) 
N,(s) = X A L l (s - S ')L 7 94( S ' + p), (14) 

s' ~ 

where SpJs' + fx) is a plaquette variable satisfying d'^S^Js) = J 7 (s) and the coordinate 
displacement fx is due to the interaction between dual variables. 

However the expectation values of the above three operators are not completely equiv- 
alent. When we consider infrared effective abelian theories, it is natural that the static 
potential between electric charges becomes Coulombic in the deconfinement phase. The 't 
Hooft operator in the dual abelian Higgs model or the Wilson loop in the generalized Villain 
form reproduce this behavior. However, it is stressed that all three operators give the same 
area law, since the differences give only Coulombic or Yukawa potentials. Since we are in- 
terested in the string tension, let us consider the operator fljjD from now on. See Appendix 
B for details. 
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D. Analytic blockspin transformation 

We construct a block spin transformation ([]) of monopole currents. [] Integrating out 
the monopole current variable on the fine lattice we arrive at an effective action and the 



loop operator for the static potential on the coarse lattice [15|. Let us start from 



(W m (C)) = £ exp \ - £ K(s)D (s - s')^') + 2m £ N^k^s) 



fcfl(s) = -oo ^ s,s',[i 



X 



II 6(K,(s^)-B k ^))/Z[k]. (15) 



The cutoff effect of the operator ( |T5"1) is 0(a) by definition. This 5-function renormalization 
group transformation can be done analytically. Taking the continuum limit a —>■ 0, n —>■ oo 
(with b = na is fixed) finally, we obtain the expectation value of the operator on the coarse 
lattice with spacing b = na 



(W m (C)) =exp{-vr 2 [°° d 4 xd A y £ jV„ (x)D \x - y)N^y) 

{ J -°° [M 

+vr 2 o 8 £ Bp(bsW)D^(bsW -bsW)B v (bsW) 



oo ( 

x y. ex P \-s\KM n) )\ 

b 3 Ku(bs) = -ao I 



+2mb 8 B,(bs^)D^(bs^ - bsW)K v (bsW) / £ Z[K, 0], (16) 

S (»), S W J ' i> 3 K M (i,s) = -oo 



where 



fUfcW) = lim a 8 2 II^(6s (n) - a«) { V - £>o~V - as')^(V), (17) 



Y /n— 1 



IL^(6s n - as) = -^5 (nas^ + (n - l)a - as M ) x J| ( £ <J (nasj n) + /a - as;) J . (18) 

"" ' ' i&n) \/=0 / 



SlK^s^)] denotes the effective action defined on the coarse lattice: 

S[K li ( S W)] = b» £ Y K »( bsin] ) D pM n] -bs inY )K u (bs( nY )- (19) 

s (n) .(„)' Ml" 



1 Note that the current K^(s^) on the coarser lattice with a lattice distance b = na satisfies the 
current conservation = by definition. 
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Since we take the continuum limit analytically, the operator ([16]) does not have no cutoff 
effect. 

The momentum representation of D^ibs^ 1 ' — bs^ ) takes the form 



1 PyPv r i(p u -y v )/2 



2^2 



(20) 



where A'^ V F 1 (p) is the gauge-fixed inverse of the following operator 



4 oo \ r 

ohii e ^(p+ztto 

^i=l li=—oo j \ 

(nti^ x2 



(p + 27r/) M (p + 27r/) t 
E,.(p + 2vr/)2 



(p + 27rQ M (p + 27rZ) l/ 
ni(p + 27r0f 



x- 



(21) 



The explicit form of D /J/l/ (p) is written in Ref. [15|| . Performing the BKT transformation 
explained in Appendix B on the coarse lattice, we can get the loop operator for the static 
potential in the framework of the string model: 



1 00 f 

(W m (C)} = (W m (C)) d x - £ exp - 7T 2 ^ ^(s)^^^- 1 ^ - Sl )Al 2 ( Sl - s')^(s') 

f/.ii/(s) = — oo ^ s,s' 

s Kf]W=° %% 



s,s> ) 



(22) 



{W m {C)) c i is defined by 

(W m (C)) d = exp I - vr 2 J°£xd 4 yJ2N^x)D^(x - y)N lx (y)\. (23) 



E. The on-axis case 



In the above calculation, we have introduced the source term corresponding to the loop 
operator for the static potential on the fine a lattice and have constructed the operator on 
the coarse b lattice by making the blockspin transformation. To check the validity of our 
analysis, it is to be emphasized that the same string tension for the flat on-axis Wilson loop 
can be obtained for I, T — > oo when we consider a naive Wilson loop operator on the coarse b 



lattice instead of that on the fine lattice (\13j). When we consider only quadratic interactions 
for the monopole action, we get the classical string tension from the large flat Wilson loop 
as follows fl5| : 



d 2 p 



sin 2 ^D-\k h fc 3 , 0, 0; i) + sin 2 y/T 1 ^, k 2 , 0, 0; 2) 



(24) 



-* (2tt) 2 

where D denotes the coupling of the monopole action determined numerically on the coarse 

oo, we can easily show that o~l agrees exactly with the 
r5|. Therefore, our analysis is natural as long as the 



b lattice. For / — >• oo and T — 
string tension derived later from 
quadratic monopole action is a good approximation in the IR region of SU(2) Gluodynamics. 
Note that we can show both quantum fluctuation parts also coincide. 



9 



IV. ANALYTICAL RESULTS OF SU(2) GLUODYNAMICS 



A. parameter fitting 

As shown already, the (numerically obtained) effective monopole action for SU(2) Gluo- 
dynamics in the IR region is well dominated by quadratic interactions. Hence we regard the 



renormalization flow obtained in Subsection HID as a projection of RT to the quadratic- 



interaction plane as written in Figure |1 We adopt the perfect operator discussed in the 
previous section as the correct one on the coarse b lattice in the low-energy SU(2) Gluody- 
namics. In order to know the explicit form of the operator, we need first to fix D (s — s'). 
This can be done by comparing D^ibs^ — bs^') with the set of numerically obtained 
coupling constants of the monopole action {Gi(b)} in Sec.|J. 

We assume D (s—s') in the monopole action (|T2j) to take aS SjS '+/3Aj i 1 (s— s')+ , jAl(s— s'), 
where a, f3 and 7 are free parameters. We can consider more general quadratic interactions, 
but as we see later, this choice is sufficient to derive the IR region of SU(2) Gluodynamics. 

The inverse operator of D (p) = a + (3/p 2 + 7£> 2 takes the form 

D o\p) = 4-^- 2 --^-X (25) 
\p A + mf p z + m 2 j 

where the new parameters k, m\ and m 2 satisfy n(m\ — m\) = 7 -1 , mf + m\ — a/7, m\m\ = 

Substituting Eq.(p5l) into Eq.(|2l|) and performing a First Fourier transform(FFT) on the 
16 4 lattice for the several input values k, mi and ni2 we calculate D^ v (p). Then one can 
obtain distance dependence of the D^ips^ — bs ( - n ' ). By matching the distance dependence 
of the D^ips^ — bs^ ) with numerical ones, one can fit the free parameters k, mi and 
W2. We find that the ratio mi/m2 is around 10 4 , but mi and m 2 can not be fixed well 
separately. Their optimal values for b = 2.1,2.9 and 3.8 are given in Table |, where we fix 
mi = 1.0 x 10 4 and m 2 = 12 for all b. The coupling constants with the optimal values are 
illustrated in Figure [|. Note that, in this figure, the lattice monopole action obtained from 
the continuum by analytical blocking also show the direction asymmetry. 



B. The string tension 



Let us evaluate the string tension using the perfect operator Q2"i2]). The plaquette variable 
S^p in Eq.(|14|) for the static potential V(Ib, 0, 0) is expressed by 

S J af s(z) = S al 8^5(z 2 )8(z 3 )e(z 1 )9(Ib - zi)6(z 4 )6(Tb - z 4 ). (26) 



In the Subsection [II B| we have seen that the monopole action on the dual lattice is in the 
weak coupling region for large b. Then the string model on the original lattice is in the 
strong coupling region. Therefore, we evaluate Eq.(^) by the strong coupling expansion. 
The method can be shown diagrammatically in Figure ^|. 
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1. The classical part 



As explicitly evaluated in Ref. [IS], the classical part of the string tension coming from 
Eq. (H) is 

cr d = —\vL — . (27) 
2 m 2 



(?ci/vphys using the optimal values k, m\ and m% are given in Table [TI], where a p h ys is the 



physical string tension. The scaling of y<J c i/a phys for physical length b seems good, although 
its absolute value is larger than 1. The difference will be analysed later in Section |V|. 



2. Quantum fluctuations 



The next to leading quantum fluctuation term comes from the second part of Eq. (p2|) . 
It corresponds to the second figure in Figure. |] and becomes [|15j 

^/ = -^e-^ 2 , (28) 

where 11(0) is the self coupling constant of the string action (P2"|). The total string tension 
is the sum a to t = cr c i + e q f- 

The quantum corrections for the string tension are given in Table |TTT| . We see they are 
negligibly small in IR region of SU(2) Gluodynamics. We can evaluate physical quantities 
using the classical part alone in the strong coupling expansion of the string model. Therefore, 
the strong coupling expansion works good and it is found that the classical string tension 
in the string model is near to the one in quantum SU (2) Gluodynamics. 



3. The on-axis case 



We evaluate next the string tension using Eq. (|24]) , where D~ 1 (k) are determined from 
the numerical data of coupling constants. By using a First Fourier transform(FFT) on the 
32 2 lattice, we perform the integration with respect to the momentum in (0). The results 
are given m Table [V|. We find that these are almost the same as those in Table U The 
validity of our analysis in Section [TT1| is confirmed. 



4- On the continuum rotational invariance 



We here comment on the continuum rotational invariance of the quark-antiquark static 
potential. For the sake of convenience we place a pair of static quark and anti-quark at 
the point (0,0,0) and (xi,X2,0) on a three dimensional timeslice, respectively. Both of 
the coordinates x\ and x 2 denote the sites sitting on the b = na lattice. Therefore the 
potential becomes dependent only on two coordinates, V = V(xi,X2). In the framework of 
our analysis [IS], the static potentials V(lb,0) and V(Ib,Ib) can be written as 
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V(lb,0) = — h.^, (29) 
K(J6,/6) = ^^ ln^. (30) 

The potentials from the classical part take only the linear form and the rotational invari- 
ance is recovered completely even for the nearest 1 = 1 sites. The recovery of the continuum 
rotational invariance of the static potential is naturally expected also for the quantum fluc- 
tuation, since we have introduced the source term corresponding to the Wilson loop on the 
fine a lattice and we have taken the continuum limit a — > 0. 



C. The glueball mass 

The mass spectrum in SU(2) Gluodynamics can be obtained by computing the cor- 
relation functions of gauge invariant local operators or Wilson loops, and looking for 
the particle poles. For examples, one can consider a two point function of an operator 
0(t) = J2xTr(F 2 )(x,t). For large time t it is expanded as 

(O(t)O(0)) ~ ^c J ,exp(-M i t), (31) 

i 

where M, is a glueball mass. 

We consider here the following U(l) singlet and Weyl invariant operator 

(t) = L-*' 2 £ Re (tt 12 + ^ 23 + *3i) (x, t) (32) 

x 

on the a- lattice at timeslice t. Here ^ij(x, t) is an na x na abelian Wilson loop and L stands 
for the linear size of the lattice. One can check easily that this operator carries ++ quantum 



number [20|. The connected two point correlation function of \I/ is given by 



(*(t).*(o)> c = (*(t).*(o))-(*(t)><*(o)> 



A £ [{<*«(*, t) ■ * 12 (y, 0)> + (* 12 (f , t) ■ ** 12 (y, 0)> 

x,y 

-2<* 12 (x, t)) 2 } + 2 {(*3i(x, *) • *i 2 (y, 0)) + (*3l(3 *) • *UV, 0)) 
-2(* 3 i(x,*))-<*i 2 (y,0))}] (33) 



Then we evaluate each expectation value in Q53D by using the string model just as done in 
the case of the calculations of the string tension. It turns out that the quantum correction 
is negligibly small and the classical part of the expectation value of the operator 0j {<D\ = 
*i 2 (x,t) • *i 2 (y,0), 2 = *i 2 (f,t) ■ ^ 2 (y,0), O s = ^i 2 (x,t), C 4 = *3i(^,*) • *i 2 (j/,0) and 
05 = \l>3i(x, t) • ^* 2 (y, 0) ) in the string representation becomes 

(0^ = exp | - 7T 2 J2^d 4 yJ2N,(x)D^(x - y)N,(y)^ (34) 

corresponding to Eq . (|23|) . 
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The plaquette variable S a p in Eq.(p4|) for (Oi)^ is expressed by 

S^(z)=S^(z)+S^(z), (35) 

Sapi 2 ) — $<xi8pzQ{azi - ayi)6{ayi + b - az 1 )Q[az 2 - ay 2 )6(ay 2 + b- az 2 ) 

5(az 3 - ay 3 )5(az 4 ), (36) 

(2) 

^a^{ z ) = ^ai^/32#(o.zi — ax\)9\ax\ + b — az\)d{az 2 — ax 2 )9(ax 2 + b — az 2 ) 

5(az 3 — ax 3 )5{azi — at). (37) 

This operator is shown diagrammatically in Figure [7]. 

Substituting this into Eq . fl34|) , one finds in the momentum representation 

= exp {-167T 2 / 7^i(e w + e ^+ i P4 t )( e -«+ e -if-s-i P ^ 



(27T) 



Since we study large b behaviors, we use the following formula 



(38) 



Then we obtain 

dp 3 dp 



/sin ab\ 2 , . . _ 

lim = nb5(a). (39) 



(Ox)* ~exp|-2fi:7r 2 6 2 y 



(27r) 2 \ p| + p§ + m\ P4 + P3 + raj 



/^P3 g ^P3 ^ J 
^T-g — cosp 3 (x 3 - y 3 )\ , (40) 

where _E P3 = yp 2 + m 2 . Since m x >> m 2 , we have neglected the term proportional to 
e -y/p*+m\t in Eq.(|g). 

Next the plaquette variable S a p in Eq.(|4]) for the (02) m * s expressed by 

S a ,(z) = -S^(z)+sS(z). (41) 

The same calculation yields 

dp 3 dp4 ( 1 1 



(0 2 )^^exp|-2K7r 2 6 2 | 



(27r) 2 \ p| + p| + m\ pl + p 3 + mj 



dp 3 e 



E P3 t 



i "^f ^ — cos^ 3 (^3 - j - (42) 

The plaquette variable S a p in Eq.(|34"D for the (O^^ is S^l(z) in Eq. (|37j) and the result 
becomes 

<°>* - exp {- w » 2 I w (rfT^i - } • (43) 
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For the opertor O4, a naive choice of S a p in Figure |8| does not contribute. But when 
S a p is chosen as in Figure ^§ the classical part (|34|) become non-zero and it is the leading 
contribution. The plaquette variable S a p in this case is expressed by 

S aP (z)=S^(z)+S^(z), (44) 

^ap( z ) = ^ai^/32^( ( 3'-2i — axi)6(axi + b — CLZi)6(az 2 — ax 2 )6(ax 2 + b — az 2 ) 
(S(az 3 — ax 3 ) — S(az 3 — ax 3 — b)) 5{az 4 — at) 
+8 a2 5^ 3 6{az 2 — ax 2 )9(ax 2 + b — az 2 )Q[az 3 — ax 3 )9(ax 3 + b — az 3 ) 
(— 5(azi — ax\) + 5(azi — ax\ — b)) 5(az 4: — at) 
—8ai8i339(azi — axi)9(axi + b — az\)6{az 3 — ax 3 )9(ax 3 + b — az 3 ) 
5(az 2 — ax 2 — b)5(az4 — at). (45) 

This leads us to 



(0 4 )t ~ exp \-nn 2 b 2 [ ^ {5 - 4cos P3 x 3 } 



+ m 2 pl+pl + m^ 



2 1.2 

-K7C b 



rjfi p ^P3^ 1 

1 * / 1 — e - i P3 x 3^ e -iP3(^3-y.i) _|_ H _ e iP3X3^ e iP3(x3-V3) 



2tt E P3 2 



(46) 



(V(t) ■ tt(0)) c ~ £ (e~ 2A+B + e- 2A - B - 2e~ 2A + 2 e - A ' +B ' + 2e^- B ' - Ae~ A ') , (47) 



Finally, we get 

Iv 

where we define 

A _ 2,2 f dp 3 dp 4 / 1 1 \ 2 

A = KTC b I — — 2 - 2 = <7 d • 6 , 

(27r) J \P4 + P3 + m 2 PI + pi + mf J 

dp 3 e~ L "-' 



r> - 2,2 f d P3 e VZ , \ 

B = kit b / — cosp 3 (x 3 - y 3 ) 

J Ztt rjp 3 



dp 3 dp4 ( 1 1 



A' = Kix 2 b 2 [ .^f {5 - 4cosp 3 x 3 } , ,, ., ,, ,, ,, ,, , . 

J (2tt) 2 1 \pl+P 2 3 + m 2 2 pl+p 2 3 + ml)' 

/fjf) Q p ^P3^ 1 r -1 
Fi ._(]__ g-*P3^3^ e -«P3(^3-2/3) _|_ ^ _ e iP3X3^ e iP3(x3-V3) 

27r i?,,, 2 L J 



Since -B and 5' contains e £p 3* ; it become very small when t » 1. Then one can expand 
the exponential and obtain finally for t » 1 



00 27T p 3 + ml 



<*(*) " *(0)>c - ^7 E (e- A 5 2 + 2e~ A 'B' 2 ) = \ L 2 e~ 2A J_ 



12 ^ _ A , f°° dp 3n n _ 1 e - 2 v / pi+^^ 



-TTFE e_A / ^{l-cospaga} „ 2 . (49) 

41/ "J^ J-00 27T P 3 + m 2 
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When t » 1, the integrand decreases repidly and the integral is well approximated by the 
suddle point value at p% = 0. Hence we get at large time t 

(tf (t) ■ *(0)} c ~ — {lV 2 ^ 2 + 2L(2L - l) e - 5fJrf fe2 } exp {-2m 2 t} , (50) 

where the second term coming from the O4 is seen to be suppressed by the factor e~ 3eTcl ' b2 , 
since a c i ■ b 2 become large for b >> 1. Other quantum corrections are also suppressed 
similarly. The lowest glueball mass M ++ is found to be M ++ = 2rri2- 

The lowest glueball mass in unit of the string tension a c i for b = 2.1, 2.9 and 3.8 are 
given in Table [IV[ This is almost consistent with the recent lattice results M(0 ++ ) / ^<J p h ys = 
3.74 ±0.12 ||. 



V. ANALYSIS 

The value of the string tension calculated analytically in the previous section is about two 
times larger than the value which is numerically determined from the monopole contribution 
to the abelian Wilson loop and is used here to fix the physical scale. 

Let us analyse the origin of the difference in details. The method and the assumptions 
we have adopted are summarized in the following: 

1. Abelian dominance We have assumed first that after abelian projection abelian 
components alone are responsible for non-perturbative phenomena of SU(2) Gluody- 
namics in the infrared region. This assumption is based on the numerical data obtained 
in MA gauge [f§,0. Bali et al. |2^] have made a detailed test at (3 = 2.5115 and have 
confirmed the assumption of abelian dominance of the string tension is good at the 
level of 92 percent. 

2. Monopole dominance The abelian Wilson operator can be factorized into monopole 
and photon contributions. We have assumed only the monopole part is responsible 
for the string tension on the basis of the numerical analysis [§,0- The values of the 



string tension we have used are listed in Table VI . The differences are not big 



DeGrand-Toussaint (DT) definition of lattice monopole We have used DT 
monopole in the numerical evaluation done in Section [D], since we do not know an 
alternative which can be used in numerical simulations. The magnetic charge of DT 
monopole is restricted. However we have used the definition of lattice monopole with 
any integer charge which we call as natural monopole in the step of the analytic 
block spin transformation. As checked in the case of compact QED [Q, there may 
be a considerable difference between natural and DT monopoles on the fine a lattice 
for small (3 region. But the difference is expected to be decreased after block spin 
transformations, since the blocked monopole can take a wider range of charge. But we 
can not estimate the effect quantitatively in the present stage. 

Truncation and scaling In the inverse Monte-Carlo calculations and numerical block 
spin transformations, we have truncated the number of the terms in the effective 
monopole action. We have used 27 quadratic terms up to 3 lattice distances and 
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four-point and six-point self interactions, assuming short-ranged interactions are more 
dominant. Then we have performed the block spin transformation the number of steps 
of which is n = 1,2,3,4,6,8. The data seem to show roughly the scaling behavior 
expected on the renormalized trajectory. However, this step could still give rise to 
fairly large systematic errors. The scaling behavior may not be enough. Actually, the 
dominant quadratic self-coupling term G(l) at b = 2.78 (J3 = 2.0, n = 4) is around 
0.16, whereas it is around 0.09 at b = 2.87 {(3 = 2.3, n = 8). 

5. Analytic calculations Since the quadratic terms seem to be dominant in the infrared 
region, we have evaluated the physical quantities in the framework of the quadratic 
monopole action. Using the mean-field approximation, the quartic term can be ap- 
proximated by the quadratic self and the nearest-neighbor terms with an effective 
coupling 8q(b) < fc„(s) >, where q(b) is the quartic coupling constant and < fc^(s) > 
is the monopole density. The induced effective self-coupling is still by two or three 
order smaller than the original quadratic self-coupling. Hence contributions from four 
and six point interactions can be neglected safely for b > l.da h y S . Since quantum 
corrections are also very small, we have made calculations using the classical contribu- 
tions alone. The strong coupling expansion of the string model calculations is reliable. 



We show the expected coupling constants of RT for large b regions in Figure 10. The 
comparison of the three parameters a (3 7 between the expected RT and the optimal 
fit to the numerical data are plotted also in Table |V11| . 

As a results, we come to the conclusion that we have to perform Monte-Carlo simulations 
on an improved action for large b starting from the points nearer to the continuum and more 
steps of block spin transformations to reproduce the correct value of the string tension. It 
is stressed, however, that the other parts of the above procedure appear rather reliable. 



VI. CONCLUDING REMARKS 

1. In order to obtain the quantum perfect effective action of low-energy SU(2) Gluo- 
dynamics, we have performed the blockspin transformations on the dual lattice after 
abelian projection in MA gauge numerically. In the Inverse Monte-Carlo method, we 
have adopted more general form of monopole actions than the one in the previous 
study P,|TT| and have stressed the important features of the almost perfect monopole 
action. We have transformed the monopole action into that of the string model of 
hadrons by using the BKT transformation. 

2. To evaluate the physical quantities, we have considered the quadratic interaction sub- 
space for the monopole action and find the correct form of perfect operators. We 
have evaluated the physical quantities such as string tension and the glueball mass 
for SU(2) Gluodynamics using the string model of hadrons analytically. The strong 
coupling expansion works good and it turns out that the classical results in the string 
model is near to the one in quantum SU(2) Gluodynamics. Probably, it means that 
the classical string theory is a good approximation for IR Gluodynamics. 

3. To get a better fit of the string tension, we have to perform more elaborate Monte-Carlo 
simulations for large b on larger lattices. 
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APPENDIX A: 

The quadratic interactions used for the modified Swendsen method are shown in Ta- 
ble |V111| . Only the partner of the current multiplied by fc M (s) are listed. All terms in which 
the relation of the two currents is equivalent should be added to satisfy translation and 
rotation invariances. 

The higher order interactions used for the modified Swendsen method are listed in Ta- 
ble pel. 
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APPENDIX B: 



In this appendix we give various representations of the Wilson loop operator. 

The original representation 

Let us consider the generalized Villain action defined by Eq. ([!]). In this model, the 
quantum average of the Wilson loop operator is written as 

(W(C)) = |exp|i^J,( S )^( S )}| = Z[J]/Z{0], (Bl) 

[•IT +°° f ~| 

Z[J}= \{de,{s) E op -5[M]+»EW» ■ (B2) 

n s ^ nflv {s)=-oo I s,n ) 

We call this as the original representation of the Wilson loop. 

The monopole representation 

The above original representation can be transformed into the monopole representation 
exactly in the following way. Let us perform the BKT-transformation with respect to the 
integer- valued tensor n pu (s) in (|B2|) : 

+ d {il q v] {s), (B3) 



d[^m vp] (s) = -e^pxkxis - p-P- p), (B4) 

where m pv (s) and q^s) are rank-2 tensor and vector fields on the original lattice respectively. 
The vector field k p (s) which can be interpreted as a monopole current on the dual lattice 
obeys conservation law d^k^s) = by definition. 
Using the Hodge-de-Rahm decomposition we write 

d^B v] (s) + 27rrv(s) = d { ,9^ c \s) + 2tt E & p &Z l (s - s') l -e ppvX k x {s' -p-fr-V) (B5) 

s' " 

ef c \s) = 6,(3) + 2ttJ2 A I 1 ( s - sXm^OO + q,(s). (B6) 

a' 

Substituting Eq. ([B5|) in Eq. (|B~2|) and integrating out the noncompact field dj™' c '(s) we get 

+00 +00 ( / 1 \ 

Z[J]= E E exp -tt 2 E W )(s-s')J,(s') 

- E ^(s)^o(s-sO^ , )-27rzE^(s)A£ 1 ( S - S / )^^( s ')}- (B7) 

s,s';fj, s,s' ) 

It is convenient to define the plaquette variable S^(s) from the abelian integer-charged 
electric current J 7 (s) by the following relation 

& p S J M {s) = J 7 ( a ). (B8) 
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By this definition, S^s) can be interpreted as the surface which is spanned on the contour 
Jry(s). The third term on the exponential in Eq.([B7]) can be rewritten as follows: 



Y,Si P {s)d p /\ L l d' u m, u {s) 

S,fX 

E^^^K^I-E^^^^l'^hW' (B9) 



S,fl S,/J, 

When use is made of Eq. (|B4|) , we have 



*ENp{a)kp{8), (BIO) 



where N^s) is defined by flU]) . 

The summation with respect to the integer field m pv (s) is trivial since 

exp{27u x integer} = 1. Therefore, the expectation value of the Wilson loop operator 
in the monopole representation becomes: 

(w(c)) = (w(c)) m • exp j-^ 2 E E U*) (a^) (* - J , (BH) 

where (W(C)) m is written as 

(W(C)) m = /exp|27riX;MaWa)l\ = S[J]/S[0], (B12) 



^1 = E II <W)>° ex P - 5 ^ + 2m E U*Ws) . (bis) 

fe M ( S ) = -00 V S / { 8,/li J 

The monopole action S[k] is shown in (0). 

Ao£e t/iat t/ie difference between (W(C)) m and (W(C)) is only an electric- electric cur- 
rent J — J interaction which comes from the exchange of regular photons and has no line 
singularity leading to a linear potential. Hence the term of the area law of both operators 
are completely the same. So concerning the low-energy physics of QCD, such a term is not 
so important. We, therefore, neglect J — J interactions and consider (W(C)) m to evaluate 
the static potential. The analysis in [IS| leads to Eq.( 



The dual representation 

As is written in the theory described by the monopole action is given in the 
particle representation. It can be expressed in the field representation as a field theory. 
This is a dual abelian Higgs model. We show here the above monopole representation is 
equivalent to the lattice form of the modified London limit of the dual abelian Higgs model. 
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Introducing an auxiliary dual field 0(s) for the constraint of the monopole current 
<5a'fc M (s),o and a dual vector field C M (s), Eq. ([B13| ) is rewritten as 



- «V„M 



+i^ l (C IM (s) + d'^(s)-27rN IM (s))k li (s)- £ ^(a)^^ - a') k,(s')\. 

(B14) 

Inserting the unity 1 = VF5 (F^s) — k^s)) to Eq. ( |B14| ) and performing the Gaussian 
integration with respect to the F^s) field, we have 

/OO °° f 1 / 

\{dC,{s) n#) E exp -- 53 (^a](5)-27r^( S 

"°° s;n J 77 s;it 1 ( s )=-oo \ 4 ^ s;fi>u 



-\ E ipM + W s ) + 2 ^( s )) ^i - ^ - *o (c^o + wo + 

(B15) 

where we have used also the Poisson summation formula 

00 00 ( \ 

£ 5 (F,(s) - k,(s)) = £ ex P j^E^WU*) • ( B16 ) 

Therefore the expectation value of the Wilson loop operator in the dual representation 
becomes 

(W(C)) = exp j-vr 2 £ ^)(«-^Jw)), (B17) 

where -ff(C) is a 't Hooft loop operator defined by Eq.(ll). We see 
(H(C)) = Z[Sj]/Z[% 

/oo pn 00 [ 1 . 

Y[dC»(s) lld^s) exp - — £ (d^C u] (s)-2nSl u (s; 

"4 E + W*) + ^(s - 5') (C M (s') + ^(s') + 2vr/,(,'; 

(B18) 

Eq.fllOD is the lattice form of the modified London limit of the dual abelian Higgs model. 
C M (s) and 4>(s) can be interpreted as a dual abelian gauge field and the phase variable of 
the dual Higgs field, respectively. Note that the integer-valued field Z M (s) appears due to the 
compactness of the theory. 
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The string representation 

We show here the string representation is obtained from the monopole representation. 
Introducing an auxiliary field <f>(s) for the constraint of the monopole current 8qi kp(8),o an d 
inserting the unity 1 = /f^ VF5 (F^s) — fc^(s)) into Eq. (|B13| ), it is rewritten as 



(s - s') J M (s' 



/ + OO /*7T °° ! 



iY F M ( d 'A s ) + 2 ^( s ) + 2 ^(*: 



-1 

S.fl 



(B19) 



Here we also have used the Poisson summation formula for the integer valued vector field 



Ms). 



Now we perform the BKT transformation with respect to the integer valued vector field 

= s^s) + d,j,r(s), d^s u] (s) = ^e^ af3 a af3 (s -&-$), (B20) 

where r(s) is a scalar field defined on the dual lattice and the string field a a p(s) defined on 
the original lattice obeys the conservation law d 1 a^s) = by definition. This means the 
variables <J aj 3(s) form a closed surface on the original lattice. 

Integrating out all fields except for the string field (Tapis) , we obtain the following string 
representation defined on the original lattice: 

Z[J] = exp I -vr 2 £ J^s) (s - s')J»(s') - vr 2 £ N,(s) (i-) (s - s')N„(s>) 



X 



E (n^M,M,o) expi-vr 2 ^ a^(s) ( -^—) (s - sO<v( s ') 



(B21) 



22 



G1 vsb 



FIGURES 



G2,G3 vs b 



2.0 



1.0 



■42* monopole (24 4 lat.) 
► 3 3 monopole (24 4 lat.) 
■ 4 3 monopole (48 4 lat.) 

♦ 6 3 monopole (48 4 lat.) 

• 8 3 monopole (48 4 lat.) 



0.0 



0.0 1.0 



2.0 



3.0 



4.0 




■*2 3 monopole (24 4 lat.) 
► 3 3 monopole (24 4 lat.) 
■ 4 3 monopole (48 4 lat.) 

♦ 6 3 monopole (48 4 lat.) 

• 8 3 monopole (48 4 lat.) 
<l 2 3 monopole (24 4 lat.) 
^3 3 monopole (24 4 lat.) 
□ 4 3 monopole (48 4 lat.) 
O 6 3 monopole (48 4 lat.) 
O 8 3 monopole (48 4 lat.) 



0.00 



0.0 



1.0 



2.0 



3.0 



4.0 



G10 vs b 



0.15 



0.10 



0.05 




-4 2 s monopole (24 4 lat.) 
► 3 3 monopole (24 4 lat.) 
■ 4 3 monopole (48 4 lat.) 

♦ 6 3 monopole (48 4 lat.) 

• 8 3 monopole (48 4 lat.) 



0.00 



0.0 



1.0 



2.0 



3.0 



4.0 



4pt. vs b 



0.10 



0.08 - 



0.06 - 



0.04 



0.02 



•*2 3 monopole (24 4 lat.) 
► 3 3 monopole (24 4 lat.) 
■ 4 3 monopole (48 4 lat.) 

♦ 6 3 monopole (48 4 lat.) 

• 8 3 monopole (48 4 lat.) 



0.00 



0.0 



1.0 2.0 



3.0 



4.0 



FIG. 1. The couplings of quadratic interaction term and 4-point interaction term versus 
physical length b. 
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FIG. 2. The renormalization flow on the projected plane. 
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Monopole action R-dependence 

b=2.14 



O (3=2.4 , 48 4 lattice , 8 3 monopole 

1.0e+00 r 

1.0e-0lT .G2 

Q^G3 

CD 

1.0e-02 ~ O 

CP O o 



1.0e-03 



1.0e-04 



G1 G2 G3 



012345678 
distance R 2 

FIG. 3. The distance dependence of the couplings of quadratic interaction terms at b = 2.14. 
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Couplings vs distances 

48 4 lattice, 8 3 monopole 
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FIG. 5. The coupling constants with the optimal values k, mi and ni2 for b = 2.1, 2.9 and 3.8 
from the comparison with numerical data. 




FIG. 6. The strong coupling expansion of the Wilson loop calculation. 
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Couplings vs distances 

48 4 lattice, 8 3 monopole 
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FIG. 10. The expected coupling constants of RT (star) versus numerical data. 
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TABLES 



TABLE I. The optimal values k, mi and m<i for b = 2.1, 2.9 and 3.8 from the Inverse 
Monte-Carlo method 

b 2.1 2.9 3.8 

~k L76 3T2 483 

mi 1.0 xlO 4 1.0 xlO 4 1.0 xlO 4 

m 2 12.0 12.0 12.0 



TABLE II. Ja d /a phys for b = 2.1, 2.9 and 3i 



6 2.1 2.9 3.8 

L64 L56 L45 

Vphys 



TABLE III. The leading quantum correction for b = 2.1, 2.9 and 3.8. 

b 2.1 2.9 3.8 

£ e -4n(o)b* 1.26 x 10- 5 1.40 x lO" 9 1.65 x lO" 14 



TABLE IV. M 0+ + /y/dci for 5 = 2.1, 2.9 and 3.8. 



2.1 

5.56 



2.9 
4.18 



3.8 
3.36 



TABLE V. Ja L /a phys for b = 2.1, 2.9 and 3.8 from Eq.(||). 



^phys 



2.1 2.9 3.8 

1.73 L59 L39 



TABLE VI. String tensions from non-abelian (<tj) and monopole (a m ) Wilson loops. 



2.20 0.4690(100) 0.4804(52) 

2.30 0.3690(30) 0.3589(36) 

2.40 0.2660(20) 0.2678(82) 

2.50 0.1905(8) 0.1851(32) 

2.60 0.1360(40) 0.1346(39) 

2.70 0.1015(10) 0.1016(21) 
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TABLE VII. The comparison of the three parameters a (3 7 between the expected RT and the 
optimal fit to the numerical data. 








1 
Z.L 


a 

2.9 




Q Q 
3.8 


a 




0.565 


0.321 




0.207 


P 




6.78 


3.85 




2.49 


7 




5.65 x 10" 5 


3.21 x 10 


-5 


2.07 x 10" 5 


OtRT 




1.52 


0.780 




0.435 


Prt 




6.78 


3.85 




2.49 


1RT 




4.09 x 10~ 4 


1.90 x 10 


-4 


9.15 x 10~ 5 


TABLE VIII. The quadratic interactions used for the modified Swendsen method. 


coupling {Gi 


} distance 


type 


coupling {Gi} 


distance 


type 


Gi 


(0,0,0,0; 




G15 


(2,1,1,0) 


k^(s + 2p + v + p) 


G 2 


(1,0,0,0; 


k^(s + /t) 


G\6 


(1,2,1,0) 


k^(s + p + 2v + p) 


G 3 


(0,1,0,0; 


k^(s + v) 


Gn 


(0,2,1,1) 


k^(s + 2z> + p + a) 


G A 


(1,1,0,0; 


k^(s + p + C>) 


Gis 


(2,1,1,1) 


k^(s + 2fi + v + p + a) 


G 5 


(0,1,1,0; 


k^s + v + p) 


G±g 


(1,2,1,1) 


k^(s + p + 2z> + p + a) 


G 6 


(2,0,0,0; 


fc M (s + 2/x) 


G20 


(2,2,0,0) 


k^{s + 2p + 2i>) 


G 7 


(0,2,0,0; 


k fl {s + 2i>) 


G21 


(0,2,2,0) 


k^(s + 2u + 2p) 


G 8 


(l,l,l,i; 


k/j,(s + p + i> + p + a) 


G22 


(3,0,0,0) 


kfj,(s + 3p) 


G 9 


(1,1,1,0; 


k^(s + /t + v + p) 


G23 


(0,3,0,0) 


kfj,(s + 30) 


Gw 


(0,l,l,i; 


k/j,(s + z> + p + a) 


G24 


(2,2,1,0) 


k fl {s + 2p + 2v + p) 


G11 


(2,1,0,0; 


k^(s + 2p + u) 


G25 


(1,2,2,0) 


k^{s + p + 2v + 2p) 


G12 


(1,2,0,0; 


k^{s + p + 2u) 


G26 


(0,2,2,1) 


k^(s + 2z> + 2/5 + a) 


G13 


(0,2,1,0; 


k^{s + 2u + p) 


G27 


(2,1,1,0) 


k p (s + 2p + 2v + p) 


Gu 


(2,1,0,0; 


k v (s + 2p + u) 








TABLE IX. The higher order interactions used for the modified Swendsen method. 


coupling 




distance 






type 


4-point 




(0,0,0,0) 








6-point 




(0,0,0,0) 






e s (e*=_ 4 W) 3 
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